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Abstract 

The energy eigenvalues of the class of non-Hermitian 'PT-symmetric Hamiltonians H = + 
x^{ixy (e > 0) are real, positive, and discrete. The behavior of these eigenvalues has been studied 
perturbatively for small e. However, until now no other features of H have been examined per- 
turbatively. In this paper the small-e expansion of the C operator and the equivalent isospectral 
Dirac-Hermitian Hamiltonian h are derived. 
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I. INTRODUCTION 



For non- Hermit ian "PT-symmetric Hamiltonians it has been established that the physical 
requirements of spectral positivity and unitarity can be met even though the Hamiltonian 
is not Hermitian in the Dirac sense. (A Hamiltonian H is Hermitian in the Dirac sense if it 
satisfies H = W , where the Dirac adjoint-symbol f indicates combined complex conjugation 
and matrix transposition.) Many "PT-symmetric model Hamiltonians have been studied 
but the first non-Hermitian PT-symmetric Hamiltonian for which spectral positivity and 
unitarity were verified is 

H = p^ + x\ixy (e>0). (1) 

It was shown in 1998 that the spectrum of the class of Hamiltonians ([T]) was positive 
and discrete ^ and it was conjectured that spectral positivity was a consequence of the 
invariance of H under the combination of the space-reflection operator P and the time- 
reversal operator T. Three years later, a rigorous proof of spectral positivity was given 
jsl. Then, in 2002 it was demonstrated that the time-evolution operator U = e~*^* for the 
Hamiltonian ([1]) is unitary [J] . In Ref . [3] it was shown that if the PT symmetry of a non- 
Hermitian Hamiltonian is unbroken, then it is possible to construct a new operator called C 
that commutes with the Hamiltonian H. The Hilbert-space inner product with respect to 
the CPT adjoint has a positive norm. Also, the operator U is unitary with respect to the 
CPT adjoint. Thus, Dirac Hermiticity of the Hamiltonian is not a necessary requirement of 
a quantum theory and unbroken PT symmetry is sufficient to guarantee that the spectrum 
of H is real and positive and that the time evolution is unitary. 

In subsequent papers the C operators for various quantum-mechanical and field-theoretic 
models were calculated SBBHi, mostly by using conventional perturbative methods. 
It was shown that this operator has a natural form as the parity operator multiplied by an 
exponential of a Dirac Hermitian operator Q: 

C = e^P, where Q = (2) 

The operator Q vanishes in the unperturbed e — limit when the Hamiltonian becomes 
Hermitian and parity invariant. This implies that the C operator can be interpreted as the 
complex extension of the parity operator P. It was proved by Mostafazadeh that the Q 
operator can be used to transform the non-Dirac-Hermitian Hamiltonian H to a spectrally 
equivalent Dirac- Hermitian Hamiltonian h (lo| : 



h = e^l^He-^'\ (3) 
This similarity transformation was used by Geyer et al. to convert Hermitian Hamiltonians 



to non-Hermitian Hamiltonians 11 



Originally, Bender and Boettcher introduced the Hamiltonian ([T]) to examine the conjec- 
ture by Bessis and Zinn- Justin that the spectrum of the Hamiltonian H = + ix^ might 
be real. Bender and Boettcher speculated that if this conjecture were true, then the reality 
of the spectrum might be due to the obvious symmetry of this Hamiltonian under combined 
P and T refiection. To study this conjecture Bender and Boettcher considered the Hamil- 
tonian in ([T]) because (i) this Hamiltonian is PT symmetric for all real e, and (ii) it could 
then be studied perturbatively for small e by using the methods of the 6 expansion, which 
had been developed much earlier by Bender et al p^]. (One recovers the Bessis- Zinn- Justin 
model Hamiltonian by setting e = 1.) The discovery that order- by-order in powers of e the 



2 



eigenvalues of if in ([T]) are all real led to much subsequent numerical and analytical work 
on this model. 

Surprisingly, the methods of the 6 expansion, for which the principal idea is to introduce 
in the exponent a small perturbation parameter whose effect is to quantify how nonlinear 
a theory is, was not used in further studies of the Hamiltonian in ([T]). The objective of 
this paper is to report a new perturbative study along these lines in which the C operator 
and the equivalent Dirac-Hermitian Hamiltonian h are calculated for small e. We have 
determined the C operator to first order in powers of e, and using this result we have found 
the equivalent Hermitian Hamiltonian h to second order in e. The results can be presented 
compactly, but they reveal in dramatic fashion how complicated and nonlocal the isospectral 
Dirac Hermitian Hamiltonian h can be. 

The construction of the C operator in Ref. j3] was the key step in showing that time 
evolution for the non- Hermitian Hamiltonian ([T]) is unitary. However, the difficulty with the 
construction given in Ref. [3] is that calculating the C operator requires as input all of the 
coordinate-space eigenvectors of the Hamiltonian. This information is available in quantum 
mechanics but it is unwieldy. (In the case of quantum field theory this information is not 
available because there is no simple analog of the coordinate-space Schrodinger equation.) 

Fortunately, it ispossible obtain the C operator by solving three simple simultaneous 
algebraic equations |5|: 

C = 1, (4) 

[C, VT] = 0, (5) 

[C,H] = 0. (6) 

The first two of these equations are kinematic because they are obeyed by the C operator 
for any PT-symmetric Hamiltonian. If we seek a solution for C in the form ([2]), we find that 
these two equations imply that the operator Q{x,p) is an even function of the x operator 
and an odd function of the p operator. The third equation is dynamical because it makes 
explicit use of the Hamiltonian that defines the theory. This is the equation that we will 
solve perturbatively using the methods of the 6 expansion. 

In Sec. iniwe calculate the Q operator to first order in e for the PT-symmetric Hamiltonian 
in ([T]) and in Sec. IIIII we calculate the equivalent Hermitian Hamiltonian h to second order. 
In Sec. HV] we discuss the calculation of eigenvalues. 



II. FIRST-ORDER CALCULATION OF THE C OPERATOR 

We begin our calculation of the C operator for if in ([1]) by expanding H to second-order 
in powers of e: 

H = Ho + eHi + e^H2 + 0{e^), (7) 

where Hq = p^ + x"^ is the Hamiltonian for the harmonic oscillator. Hi = x'^log{ix), and 
H2 = |[x log(za;)]^. We then recall the representation of the C operator in ([2]) and expand 
the Q operator as a series in powers of e: 

00 

2 = E^"2„. (8) 

n=l 

This perturbation series begins at = 1 because when e = 0, the C operator reduces to the 
parity operator V. We will see that even-n terms as well as odd-n terms must be included in 
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([8]). This is a significant departure from previous perturbative results for the C operator; for 
the cubic Hamiltonian H = p'^+x'^ + eix^ ^ and the square- well Hamiltonian H = p^ + Vlx), 
where V{x) = eix/\x\ {\x\ < 1) and V{x) = oo > 1) [7], there are no even-n terms. 
If we then substitute the expansion for C, 



C 



l + eQi+e2Q2 + ^e'Q? + 0(e=^) 



and the expansion for H in into the commutator ([6]) and collect powers of e, we obtain 
a sequence of equations for the coefficients of Q. After some algebra, we find that the 
first-order equation simplifies to 

[Qi, Hq] = x^ [log(za;) — \og{—ix)] = inxlxl (9) 

and that the second-order equation becomes 

[Q2, Ho] = iiTx\x\ log |x| - [Qi, x^ log . (10) 



Although we have not yet found an analytical solution to (jlOl) . it is likely from this 
equation that Q2 is nonzero, and as we stated earlier, this is an unexpected result based on 
previous perturbative calculations of the C operator. Although ([9]) is simple looking, it is 
difficult to solve. Nevertheless, we have found an exact analytical solution. 

The solution of ([9]) relies heavily on the work of Bender and Dunne [l^]. We introduce 
the set of Weyl-ordered operators 



T 



= (m,n = 0, 1,2, ...). (11) 



The operator Tm,n is a totally-symmetric quantum-mechanical generalization of the classical 
product p'^x". Weyl-ordered operator products rely implicitly on the Heisenberg algebraic 
property that = i. We then define the generalized Weyl-ordered operator Tm,n as 

(12) 

fc=0 ^ ^ 

In terms of this definition of Weyl ordering, we assert that Qi can be expressed as 



2i - i E (2n-l)(L + l)^^"+^'-^"-^^- ^^^^ 

It is clear that Qi satisfies the kinematic constraints that it be even in x and odd in p. We 
do not claim that f|T3l) is the unique solution to ([9]). 

To show that Qi solves ([9]), we first demonstrate that the n = term in the series 
commuted with x'^ gives i7ra;|a;|; that is. 



2 L 



TlA,X^ 



mx\x\. (14) 



We then show that the /cth term in the series commuted with is exactly canceled by the 
{k + l)st term commuted with x"^. 



4 



To verify ([Hj) we note that the n = term in the series is — ^Ti^i = — | (|a;|j9 + p\x\). We 
then get 



■— {\x\p + p\x\) , = — — (|x| [p, x^~\ + [p, x^] = i7rx|x|. 



Next we show that the kth term of Qi commuted with p"^ gives 



71 

2 



{2k-l){2k + l) 



rp 2 

-'2fc+l,-2A;+l;P 



m{-lY 



22fc+2(2A;+ 1) 



2fc+2 

E 

j=0 



2A; + 2 



where the sign function sgn(x) is defined by 



p^sgn(x)a;~^V^"^^~^ 
(15) 



sgn^xj 



1 if X > 0, 
if X = 0, 
-1 if X < 0. 



The identity ^|x| = sgn(x) imphes that [|x|,p] = 2sgn(x). Also, from the identity 
^ (jxl^lxl"") = n|x|~"'"Sgn(x) + |x|"^|x|~" for all integer n, we obtain [|x|~^'^"'"^,p] = i{l 
2fc)sgn(x)x~^'^. It follows that 



T, 



2k+l-2k+liP 



1 /2k + 1 



22fc+l z i \ n 

j=0 ^ 



p'lxl 



^ 2k+l 
fc+T ^ 



22fc+l 



3=0 



2k + 1 
j 



~2k+lp2k+l-j 



^p^ [|x|-2'=+\p]p2'=+i-^' 



2fc+l 



i{l-2k) 

22k+l 2^ 
j=0 



2k + 1 
3 



pPsgn{x)x 'p 



2k 2k+l-j 



Using this result, we find that 



T2k+l-2k+l,P I = P 
2fc+l 



i{l-2k) 



22fc+l 

i(l -2fc) 

22A:+1 



. r=0 
■2fc+2 



2A; + 1 



T2k+l-2k+l,P + ^2fc+l,-2fc+l,Pj P 

2k+l 



^ /2A; + 1 
V r' - 1 



r'=l 



p''+^sgn(x)x"^V'+^^'^ 



p'''sgn(x)x 2fep2fc+2 r' _^ 



E 

s=0 



2A; + 1 



P''sgn(x)x 2fcp2A:+2 



+sgn(x)x-2V"^^ 



2fe+l 

+ E 



2A;+ 1 
s 



p sgn(xjx p 



2k'2k+2-s 



i{l-2k) 

22fc+l 



sgn(x)x-2V"^' +/"^'sgn(x)x-2" + 



2k+l 

+ E 



2A; + 1 
J-1 



2A; + 1 
j 



p^sgn(x)x p 



2k 2k+2~j 
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We then simplify this expression by using 



2k + 1\ f2k + 1\ f2k + 2 



J-1 



3 



3 



which estabhshes the result in ( 1151) . 

We now calculate the {k + l)st term of Qi commuted with x^. We first note that T„ 
can be rewritten as Ui 



T 

J. rr 



n 

-Y 

On / ^ 



3=0 



n 



Using the identity 

[p", x^] = -in {jP-~^x + xp^"^) = -2inxp^~^ + n{n - = -2inp"-~^x - n{n - l)/)""^, 



we obtain 



T T 



1 

-Y 

on / ^ 



\Y.{^\A'{-^rn{xr-^^r-^, 



x]\\x 



i=0 

— 2zm 



-2im 



i=o 



m— 1 



r 



r=0 



m — 1 



Thus, 

2fc+3,-2fc-l, 3: 

This shows that 



-2i{2k + 3) (2k + 2 

22A:+2 / > 

r=0 



TT 

2 



(2fc-l)(2fc + l) 



2fc+l,-2fc+l;P 



22fc+2(2A; + 1) 



i=o 



2A; + 2' 



p'e[x)x p 



2k 2k+2-j 



We conclude that the fcth term of Qi commuted with p^ is exactly canceled by the [k + l)th 
term commuted with x^: 



TT 

2 



-1 



(2A; + l)(2A; + 3) 



rp 2 

-'2fc+3,-2fc-l; 3; 



TT 

2 



(2A;- 1)(2A; + 1 



rp 2 

-'L2k+l,-2k+l,P 



We have thus shown that Qi in ( fT3l) solves the commutation relation ([9]). This argument 
was quite elaborate, and it is clear why we have not yet found an analytical solution to the 
commutation relation (ITU]) . However, having found Qi, we can now calculate the equivalent 
Hermitian Hamiltonian h to second order in e, as we show in the next section. 
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FIG. 1: Two plots of the potential V{x) = + ex^\og\x\, one for the unperturbed case e = 
(dashed line) and the other for e = 1 (solid line). Note that when e > 0, the potential x'^+ex^ log |a;| 
lies below for |x| < 1, but for |x| > 1 it rises faster than the parabolic potential and thus squeezes 
the energy levels upward. 



III. SECOND-ORDER CALCULATION OF THE EQUIVALENT HERMITIAN 
HAMILTONIAN 

In this section we use ([3]) to calculate the Hermitian Hamiltonian h, which is isospectral 
to H 'm. ([T]). Expanding ([3]) as a perturbation series in powers of e, we obtain 

h = e ^^e. = ^l-- + - + ...j//^l + - + - + ...j 

+6^ [Ho, Q2] + \ {Ho, Ql} + \ [H,, Q,] - ^Ho^ + H,^ + 0(6^), (16) 

where curly brackets indicate anticommutation relations. We evaluate the first- and second- 
order terms in this equation and show that they can be reduced to compact forms. 
To first order in e, we use to simplify flTBl) and get 

(ITT \ 
——x\x\ + log(z2;) j = + + ex^ log (17) 

Thus, to first order in e the potential for the equivalent Hermitian Hamiltonian /i is a minor 
correction to the potential for the harmonic oscillator. As shown in Fig. [H when e is 
positive, the potential + ex^ log |x| lies below x^ for |x| < 1, but for |x| > 1 it rises faster 
than the parabolic potential and thus squeezes the energy levels upward. 
To second order in e, the equivalent Hermitian Hamiltonian takes the form 

h = Ho + ex^ log|x| + e'^f{x,p) + O(e^), 
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which can be simphfied by using (fTOl) : 

fi^,p) = \ [Ho, Q2] + I {Ho, Ql} + I [H,, Qi] - ^Ho^ + H, 



^QlHo 



iTTx\x\ log \x\ — Qi [Qi, Hq 



+ 1 {Ho, Ql} + ^ [H^, Qi] - ^QiHoQi + H, 
= \ [Ql Ho] - [Qi, Ho] + i {ifo, Q?} - IQiHoQi 

+ii2 ^^fI log \x\. 

The significance of this formula is that we do not need Q2 to calculate h to order e^. 
The result in (1181) may be further simplified by using 

H2 — —x\x\ log |x| = -x^ [log(ix)]^ — -iiTx\x\ log |x| = -x^ (log — — 



and 



^ [Q?, ifo] - ^Qi [Qi, Ho] + i {ifo, Ql] - IQiHoQi 
= ^QlHo - ^HoQl - \qIHo + ^QiifoQi + \hoQI + ^Q?ifo - \QiHoQi 
= -}-QIHo + \q,HoQi-\hoQ\ 
= \[[Q„Ho],Q,]. 
Thus, the second-order correction can be written as 

/(x,p) = \x' (log |x|)^ - y^' + ^ [[Qi. ^0] , Qi] . 

We then calculate the double commutator [[Qi,ifo] , Qi] as follows: 
[[Qi,i7o],Qi] = ^7r[x|x|,Qi] 

- [x^, Qi] if X < 



f (2n-lK2n+l) t^^' ^2n+l -2n+l] if X > 

f Xl^O (2n-lKL+l) [^^' ^2n+l -2n+l] if X < 



^ 2^ 7^j;—j\(^j;;—^ F . ^2n+i,-2n+ij • 



_ ^^J2n-l)(2n + l) 
Recalling that [x^, Tg^+i _2n+i] = 1i{2n + l)T2„_2„+2 [isl], we get 

f(x.)-'x^(wixn^ "'v ^-'^V 

n=0 
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A major result of Ref. 13|] is that 



n=0 

where 



°° ( — 1)" / l\ 



F{x,p) = —— / ds 



2Hq Jq cosh( 



2Ho^ 



We simply take the Dirac-Hermitian conjugate of this equation, which inverts the order of 
the operators p and x. The definition of Weyl-ordered operators Tm,n guarantees that they 
are Dirac-Hermitian. Then, noting that F{x,py = F{x,p), we deduce that 

- — —T2n+i-2n-i = arctau i-p] + iF{x,p). 
^-^ in + 1 \x 

n —0 ^ ' 



n=0 

Hence, 



n=0 ^ ^ ^ ^ 



Finally, using {x,Tm.„} = 2Tm,n+i and {p,Tm,„} = 2X^+1^^ [13|, we obtain 



^-^ 2n — 1 ^-^ 2m + 1 

n=0 m=— 1 



—X^ — - Ix, Ip, \2 7-rT2m+l ,-2m-l 

2 1 f f M^ A 

— x — - < X, < p, arctan \p— I + arctan I — p 



where the curly brackets indicate anticommutators. Thus, we obtain an explicit formula for 
the equivalent Hermitian Hamiltonian h to second order in e: 



h = Ho + ex'^ \og\x\ 



l^-x^ (log + — |a;, |p, arctan (^~^ + arctan (^~pj j" j" J ^'^^^ 

Note that h is singular and nonlocal because it contains all positive powers of p and all 
negative powers of x. 

IV. CALCULATION OF ENERGY EIGENVALUES 

If we expand ([1]) as a series in powers of e, we obtain 

H = p^ + ex^ \og{ix) + ie^[log(zx)]^ + . . . , 
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e 


Numerical value 


First order 


Second order 


0.1 


1.003 097 


1.000 912 


1.003 241 


0.01 


1.00011436 


1.000 09122 


1.000114 51 


0.001 


1.000 009 355 38 


1.000 00912249 


1.000 009 355 22 



TABLE I: Comparison of the exact ground-state eigenvalue for the non-Hermitian Hamiltonian 
in ([1]) and the first- and second-order perturbative calculations of the ground-state eigenvalue as 
given in 

where log(ix) = log(|x|) + |i7r sgn(x). We can then use conventional perturbation techniques 
to calculate the ground-state energy as a series in powers of e: 

Aground state = I + ae + be^ + O (e^) , (21) 

In first-order perturbation theory the coefficient a is the expectation value of the perturb- 
ing potential ex^ \og{ix) in the unperturbed harmonic-oscillator ground state, whose wave 
function is exp(— x^/2). We find that 

a = ( - ) = i - 2 - llog2 = 0.00912249... . (22) 
4 V2/ 2 4 2 ^ ^ ^ 

The same result for a is obtained by truncating the equivalent Hermitian Hamiltonian h in 
f l20p after the first-order term. 

We calculate the coefficient b in (l?T]) by using second-order perturbation theory applied 
to H. The calculation is lengthy, and we do not discuss it here except to give the result: 

6= -^[8 + 27r2log2 + 4(7-2 + 21og2)(7 + 21og2) + 7C(3)] = 0.23289..., (23) 
128 

where 7 = 0.5772 ... is Euler's constant. 

Note that the ground-state energy of the non-Hermitian Hamiltonian H = + x'^{ixY 
is slightly higher than the ground-state energy of the harmonic oscillator. This is consistent 
with previous numerical calculations in Ref. ^ and agrees with a precise numerical calcu- 
lation of the ground-state energy for the non-Hermitian Hamiltonian ([1]) when e is small, as 
is shown in Table I. 



Because the second-order equivalent Dirac-Hermitian Hamiltonian h in (120|) is nonlocal, 
it is not clear how to use h to calculate the energy eigenvalues beyond first order. If we try 
to do so by using standard Rayleigh-Schrodinger methods, we are faced with the problem of 
calculating the expectation value of the arctangent functions in (!20|) . This is a singular and 
ill-defined calculation that requires the introduction of a regulator. The singular nature of 
this calculation can be seen from the identity ^p|0) = i\0). This identity gives the formal 
result arctan (^p) |0) = arctan(i)|0), which is infinite in the absence of a regulator. In fact, 
the singularity is even more severe than this observation suggests: If we attempt to evaluate 
the expectation value of the arctangent functions without a regulator and use the Taylor 
expansion 

00 

arctan(s) = ^(-l)'s''-+V(2A; + 1), 

A:=0 
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we can only recover the value of b in (123|) if we can interpret the sum of the divergent series 
+ l + + l + + l + + ...as| log(2). 

We conclude that, while the equivalent isospectral Hamiltonian h is formally Dirac Her- 
mitian, using h to calculate the energies presents serious difficulties. We have encountered 
here the same kind of difficulties that were discovered in Ref. [l5|, namely, that the Dirac- 
Hermitian Hamiltonian h suffers from nonlocality and consequently is hard to use in a 
calculation. In Ref. |[15i] dimensional regulation was used to eliminate divergences in Feyn- 
man diagrams. We defer to a future paper the search for a suitable regulator for the present 
problem. 
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